Zero-bias anomaly in cotunneling transport through quantum-dot spin valves 
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We predict a new zero-bias anomaly in the differential conductance through a quantum dot 
coupled to two ferromagnetic leads with antiparallel magnetization. The anomaly differs in origin 
and properties from other anomalies in transport through quantum dots, such as the Kondo effect. 
It occurs in Coulomb-blockade valleys with an unpaired dot electron. It is a consequence of the 
interplay of single- and double-barrier cotunneling processes and their effect on the spin accumulation 
in the dot. The anomaly becomes significantly modified when a magnetic field is applied. 
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Introduction. - The combination of Coulomb inter- 
action effects, which frequently are strong in nanostruc- 
tures, and spin-dependent transport in systems coupled 
to ferromagnetic leads opens a new field of research 
with qualitatively new transport properties P, Q . Spin- 
dependent transport through nonmagnetic grains may be 
influenced by the presence of spin accumulation 0| lead- 
ing to a different transmission for parallel and antiparal- 
lel orientation of the leads' magnetization, which results 
in a finite tunnel magnetoresistance 0]. In the limit of 
weak dot-lead coupling and when a dot level is in reso- 
nance with the Fermi level of the leads (linear response 
regime) , transport is dominated by sequential tunneling. 
Away from resonance, sequential tunneling is exponen- 
tially suppressed, and transport is due to higher-order 
tunneling Q . In the Coulomb-blockade valley with an 
unpaired electron occupying the dot, Kondo-assisted tun- 
neling 7] gives rise to a pronounced zero-bias anomaly in 
the differential conductance at temperatures below the 
Kondo temperature Tk. Above Tk, transport is domi- 
nated by (second-order) cotunneling, with regular zero- 
bias behavior for nonmagnetic leads. 

In this paper, we study cotunneling transport through 
a single-level and singly-occupied quantum dot attached 
to ferromagnetic leads. When source and drain electrodes 
are magnetized antiparallel to each other, we find a pro- 
nounced zero-bias anomaly that is completely unrelated 
to Kondo correlations. It is rather a consequence of 
the interplay of spin accumulation and spin relaxation 
due to spin-flip cotunneling. A finite spin accumula- 
tion on the quantum dot partially suppresses transport. 
Spin-flip cotunneling provides a channel of spin relax- 
ation and, hence, reduces the spin accumulation. As we 
show below, single-barrier spin-flip cotunneling (in the 
absence of a magnetic field) plays a role in linear re- 
sponse, \eV\ ^ kBT, but is negligible in the opposite 
limit. This gives rise to a zero-bias anomaly in the dif- 
ferential conductance. The prediction of this zero-bias 



anomaly as well as the study of its properties is the cen- 
tral issue of this article. 

Model and Method. - We consider a quantum dot with 
a single level at energy e coupled to ferromagnetic leads 
with either parallel or antiparallel magnetization direc- 
tions. The model Hamiltonian is H = Hi^ -\- -\- iJn + 
Ht. The terms Hr = X^go- ^rqaclqcrCrqa for r = L, R rep- 
resent noninteracting electrons in the left and right lead, 
where Erqa- denotes the energy of an electron with wave 
number q and spin a in lead r. The dot is modelled by 
Hb = J2A^ ± A/2)4d<T + C/d|d|d|dx, where A is the 
Zeeman energy due to an external magnetic field, U is the 
charging energy, and the -t- (— ) sign corresponds to a 
(i). Tunnehng between the dot and leads is described by 
= Y.rqA^r-clq^da + t*dlcrqa)- Fcrromagnetism of the 
leads is included via spin-dependent densities of states, 
pI ^ pI.. The degree of spin polarization in the leads is 
characterized by the factor = [pt ~ PV)l{pt + Pr)^ 
where p^ (p~) is the density of states for spin-majority 
(spin-minority) electrons. The tunnel- coupling strength 
is characterized by FJ! = 2n\tr\^p'^. Finally, we define 
F,, = (Fj + F^)/2, F'^ = F[ + F;^ , and F = Fl + Fr, and 
assume Fl = Fr = F/2. 

We consider the Coulomb-blockade valley in which the 
dot is singly occupied with either spin. The probabilities 
for occupation with spin a are • We determine the rate 
^"p'" for a cotunneling process, in which one electron 
leaves the dot to reservoir r' and one electron enters from 
r with the initial and final dot state being a and a' , 
respectively, in second-order perturbation theory. For 
a = fj', i.e., when the dot spin is not changed, and A = 0, 
the corresponding rate is 
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while we get 
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for cotunneling process in which the dot spin is flipped 
{a is the opposite spin of a). Here, f{uj — ^ir) is the Fermi 
function of reservoir r with electro-chemical potential /i^- 
The probabilities are obtained from the station- 



ary rate equation = Y.rr' PWr?"^ ^ ^1% 
with the normalization condition _P| + 
rent / is, then, given by 
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FIG. 1: (a) Differential conductance in the parallel and an- 
tiparallel configurations as a function of bias voltage for dif- 
ferent values of spin polarization p at e = —(7/2, U — SOP, 
ksT = 0.5r, and (b) for different temperatures and p — 0.5. 



Results in the absence of magnetic field. ~ We con- 
sider symmetrically polarized leads, pl — Pr = P and, 
first, ignore a Zeeman splitting, A = 0. The differen- 
tial conductance G = dl/dV = {e^ /h)g as a function of 
the bias voltage is shown in Fig.^a). For leads magne- 
tized in parallel, we find the typical parabolic behavior 
of the cotunneling conductance with increasing bias volt- 
age. This is distinctively different for the antiparallel con- 
figuration, for which a zero-bias anomaly appears. The 
bias-voltage dependence of the conductance for different 
temperatures is shown in Fig. H^b), where the width of 
the peak grows with T. It is possible to get some insight 
into the basic properties of this behavior using analytical 
results obtained in a limit in which the formulas simplify 
considerably while all the main physics of the zero-bias 
anomaly remains included. Deep inside the Coulomb- 
blockade regime, we can neglect corrections in the ratios 
A/B with A = \eVlk^T and B = \e\,e + U. We find 
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for the parallel configuration, independent of \eV\/k-QT. 
For the antiparallel configuration, we get 



[e + Uy \e\{e + U) 



in linear response, \eV\ <C k^T, and 
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for \eV\ > k^T. We see that > g^^ 
displays a tunnel magnetoresistance. 
we find a zero-bias peak, since g^ 



i.e., the system 
But, furthermore, 
whose rela- 
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five strength, characterized hy x = {g^l^ - 5mhi)/5nun> 

X — p'^ at the edges (|e| ^ e 
|e| > e+U) tox = V/(3-p^) in the middle (e ^ -U/2) 
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of the Coulomb-blockade valley. 



The zero-bias anomaly, present (for A = 0) only in the 
antiparallel configuration, has the following properties: 

(i) The crossover from g^^^x 5miii around \eV\ w 
kBT^/8/{TTW), i.e., the width of the zero-bias anomaly 
scales linearly with fceT and depends only weakly on p. 

(ii) The relative peak height x increases monotonically 
with p and when moving from the edges towards the mid- 
dle of the Coulomb-blockade valley. 

(iii) The absolute peak height g^ax ~ Smin depends non- 
monotonically on p, since it vanishes for p = and p = 1. 

(iv) At low temperature both g^^^ and g^^ increase with 
temperature, gtl^^^ jT) / gff^^ ^.jO) '= 1 + {T / Bf + 
0{T'^) with the same constant B, such that x is nearly 
independent of temperature. 

Processes responsible for the zero-bias anomaly in the 
cotunneling regime are of second order, while these lead- 
ing to the Kondo effect are of higher than second or- 
der. The zero-bias anomaly of the cotunneling current 
is therefore distinctively different from that associated 
with the Kondo effect. The latter occurs at low temper- 
ature, T< Tkj shows up in the parallel configuration as 
well 'g', '101, grows logarithmically with decreasing tem- 
perature, and reaches perfect transmission, g = 1. Its 
width at low temperature saturates at fceTk, and it has 
a different magnetic-field dependence. 

We close with the remark that the exchange field due to 
the presence of ferromagnetic leads discussed in Ref.lglllll 
does not affect transport in the case considered here. 
Spin precession does not appear since the leads are mag- 
netized collinearly. 

Mechanism of the zero-bias anomaly. - To under- 
stand the mechanism of the zero-bias anomaly we dis- 
tinguish between four different types of cotunneling pro- 
cesses. In each of them two tunneling events are in- 
volved, either through the same or through the two op- 
posite tunnel barriers. Accordingly, we refer to them 
as single-barrier [Fig. |2Ia)] and double-barrier cotunnel- 
ing [Fig. [2Ib)]. Furthermore, the two electrons involved 
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FIG. 2: Single-barrier (a) and double-barrier (b) cotunneling 
processes, and the occupation probabilities for spin-up and 
spin-down electrons in the antiparallel configuration (c). The 
parameters are kaT = O.SF, U = SOT, £ = — (7/2, and p = 0.5. 



may carry the same or opposite spin, i.e., both single- 
barrier and double-barrier events come as either spin-flip 
or non-spin-flip cotunneling. In calculation we have taken 
into account all possible cotunneling processes. Here, 
however, we discuss just the ones responsible for the 
anomaly. Double-barrier cotunneling contributes directly 
to the current, while single-barrier cotunneling preserves 
the total charge. Nevertheless, spin-flip single-barrier co- 
tunneling can influence the total current indirectly, by 
changing of the magnetic state of the dot. In the an- 
tiparallel configuration, the dot hosts a nonequilibrium 
spin accumulation m = {P^ — Pi) /I. A different occu- 
pation of up- and down-spin levels in the dot, P| ^ P|, 
appears (even for A = 0) when the spin-flip cotunneling 
rates that change the dot from T to J, and J, to t are dif- 
ferent from each other. In equilibrium, = 0, both rates 
are trivially the same and, hence, P| = ^'x- The situation 
is different at finite bias voltage and antiparallel magne- 
tized electrodes. Now, only the two spin-flip processes 
that transfer an electron from the left to the right lead 
determine the magnetic state of the dot. The one shown 
in Fig.|2Ib) changes the dot spin from i to t- Since only 
majority spins of the electrodes are involved, the corre- 
sponding rate is larger than that of the other process that 
changes the dot spin from "f to J, by using minority spins 
only. This results in a nonequilibrium spin accumulation 
TO > (P| > Pi) that increases with V [Fig. Efc)]. The 
initial state for the dominant spin-flip cotunneling pro- 
cess that contributes to the current, Fig.Efb), is [. Thus, 
the reduced probability P| decreases transport. This is 
the mechanism by which spin accumulation gives rise to 
the tunnel magnetoresistance effect, > g^^ . 

Any spin-flip process, that reduces the spin accumu- 
lation will enhance the conductance. Such a process is 
provided by single-barrier spin-flip cotunneling. The cor- 
responding rate scales with k-oT while that of double- 
barrier cotunneling is proportional to max{|ey|, A:bP}- 
This explains the zero-bias anomaly: For \eV\ < k-gT, 
single-barrier processes play a significant role, and, there- 
fore, the current increases relatively fast with applied 



bias, which yields g^^x- 1^°^ l^^l ^ k^T, on the other 
hand, the relative role of single-barrier spin-flip processes 
is negligible as compared to double-barrier cotunneling, 
and the conductance is reduced to g^^- It is, thus, 
the interplay of spin-dependent single- and double bar- 
rier cotunneling processes that gives rise to the zero-bias 
anomaly in the differential conductance. This zero-bias 
anomaly is distinctively different from experimentally- 
observed peaks in the differential conductance for non- 
magnetic systems which occur at the onset of sequen- 
tial tunneling. 

Finally, we remark that no zero-bias anomaly occurs 
in the Coulomb-blockade valleys with an even number (0 
for e > and 2 for e -f [7 < 0), as in this case the total 
dot spin is zero, and spin accumulation is absent. 

Results in the presence of magnetic field. - In the 
presence of an external magnetic field, the dot levels are 
split by a Zeeman energy A. We restrict ourselves to 
the case of an external field that is coUinear with the 
leads' magnetization directions. When the Zeeman split- 
ting A is larger than both temperature and bias voltage, 
A| ^ maxjfcBT, |ey|}, only the lower spin level is occu- 
pied, i.e., the dot is fully polarized. In this case, spin- flip 
cotunneling is completely suppressed, which leads to a 
reduction of the conductance. With the same approxi- 
mations as we used for A = 0, we obtain 
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(e-|A|/2)2^(£ + C/+|A|/2)2 
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for the parallel configuration. Here, +/~ correspond to 
the cases when the Zeeman splitting favors a dot polar- 
ization that is parallel or antiparallel to the leads, re- 
spectively. We remark that for e = —U/2 the conduc- 
tance is symmetric under reversal of the magnetic field, 
in contrast to e —U/2, were the field reversal results 
in a different conductance [Fig. ETc d)], similar to the 
spin-readout scheme proposed in Ref. Il2l For antiparal- 
lel magnetized leads, we find 
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The above expressions approximate the plateaus shown 
in Fig. 13 When the bias voltage is increased such that 
\eV\ ^ |A| 3> k^T, spin-flip processes are again possible, 
and the dot is no longer fully polarized. For antiparallel 
magnetization of the leads and fixed orientation of the ex- 
ternal field, the conductance is asymmetric under reversal 
of bias voltage, see Fig. OJb). The reason is that single- 
barrier spin-flip cotunneling favors the dot state with the 
lower energy, which is, depending on the direction of the 
voltage drop, either the right or the wrong initial state for 
the dominant double-barrier cotunneling contribution to 
the current. The magnetic-field dependence of the peak 
height reflects the fact that single-barrier spin-flip cotun- 
neling increases linearly with A. 
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FIG. 3: (color online) Differential conductance in presence of 
a Zeeman splitting with p = 0.5, ksT = Q.2V, and U = 40r 
for the symmetric (a,b) and asymmetric (c,d) Anderson model 
with parallel (a,c,d) and antiparallel (b) relative magnetiza- 
tion of the leads. In (c), the Zeeman splitting favors a dot 
polarization parallel, in (d) antiparallel to the leads. 

In the parallel configuration, a plateau evolves for 
maxj/cBT, \eV\] <^ |A|, again due to suppression of the 
spin-flip contributions to the current. In addition, find 
that a zero-bias anomaly evolves even for parallel leads, 
when the Zeeman splitting favors a dot polarization in 
the same direction, Fig.|2Ic). As the dominant transport 
channel is non-spin-flip cotunneling with majority spins, 
and for \eV\ <C |A| < fceT, thermal occupation of the 
dot's magnetic states favors occupation with a majority 
spin, transport is enhanced as compared to A = 0. 

Summary. - We predict a zero-bias anomaly in co- 
tunneling transport through quantum dots attached to 
ferromagnetic leads that are magnetized antiparallel to 
each other. This zero- bias anomaly originates from the 
interplay of single-barrier and double-barrier spin-flip co- 
tunneling processes. From an experimental point of view, 
the anomaly may be observed in quantum dots and/or 
molecules attached to ferromagnetic leads, which include 
an odd number of electrons. Such structures have been 
already realized experimentally 0, 0, . 
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